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We report an optical detector with tunable positive operator- valued measures (POVMs). The 
device is based on a combination of weak-field homodyne techniques and photon-number-resolving 
detection. The resulting POVMs can be continuously tuned from Fock-state projectors to a variety 
of phase-dependent quantum-state measurements by adjusting different system parameters such as 
local oscillator coupling, amplitude and phase, allowing thus not only detection but also preparation 
of exotic quantum states. Experimental tomographic reconstructions of classical benchmark states 
are presented as a demonstration of the detector capabilities. 
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Detecting optical fields plays a key role throughout 
physics. From the definition of the SI unit of luminos- 
ity [1 , to the characterization of quantum processes such 
as quantum logic gates [2], precise measurement of the 
electromagnetic field is central to both fundamental and 
applied physics. The usual method of measuring fields 
in optics is either by counting photons or by measuring 
the amplitude and phase of the electric field. It is usu- 
ally not possible to move continuously between these two 
measurement regimes with a single detector. For exam- 
ple, standard homodyne detectors cannot directly probe 
the particle nature of light since it is masked by the lo- 
cal oscillator and electronic noise [3 . Conversely photon- 
number-resolving (PNR) detectors possess no phase ref- 
erence and thus have no sensitivity to the wave nature of 
light [4 . 

The action of a given detector can be specified by 
a positive operator- valued measurement (POVM) set 
{n^^}, where {[3} labels the outcomes and {7} labels 
the settings, such that for each 7 the set is complete, 
^^11/3^ = I Jp. The probability of obtaining outcome 

(3 for setting 7 and input state p is = Tr[II/3^p]. In 
conventional detectors {II /37} is fixed by the intrinsic na- 
ture of the device. Typically, such POVMs can encom- 
pass either Fock-state projectors (for PNR detectors) or 
field-quadrature projectors (for homodyne devices). In 
this Letter, we report a configurable detector with a flex- 
ible POVM, able to transition smoothly from quadra- 
ture to photon-number detection. The detector is con- 
structed from a variable reflectivity {R) beam-splitter 
(BS), two PNR detectors, and an auxiliary weak coher- 
ent state acting as the local oscillator (LO). We denote 
such device as photon-number-resolving homodyne de- 
tector (PNRHD). The BS input modes, labeled a^" and 
6^", correspond to the LO and the signal (p), respec- 
tively (see Fig. 1 (a)). The output modes, labeled by 
a°^* and 6°^*, are detected by PNR detectors Da and 
Di) giving joint outcomes {[3 = (ka^kb)}^ where ka(b) la- 
bels the number of clicks registered at Da^b)- The ad- 
justable local oscillator has settings {7 = (|a|, ^)}, where 




FIG. 1: (a) Proposed scheme for the POVM configurable de- 
tector, Da and Db are PNR detectors, a is a weak coher- 
ent state; (b) and (c) POVM elements 11/3^ corresponding to 
click events (b) pi = (1, 3) and (c) p2 = (1, 1) for LO settings 
7 = (\a\ = 1.5, 6* = 0), i? = 50% and 90% efficient PNR detec- 
tors. 11/3^ projects onto a single-mode (b) Schrodinger-kitten 
state and (c) squeezed state, with high probabilities. 



a = \a\e^^ represents the LO complex amplitude. By 
tuning the LO coupling, amplitude and phase the de- 
tector POVM set can be configured to project onto a 
variety of fundamental quantum states of the radiation 
field - Fock, displaced-Fock, quadrature-squeezed, and 
Schrodinger-kittens states for example. Figures [l] (b) 
and (c) depict the Wigner representation [6 of two 
POVM elements (II/37), for two different detector out- 
comes {/3i^/32} and LO setting 7. Here {II/37} project 
onto (b) Schrodinger-kitten states and (c) quadrature 
squeezed states with high probabilities. Since the action 
of a measurement is not only to reveal some property 
of the state of a system, but also to project the system 
in a state commensurate with that information, all mea- 
surement devices may in principle be used as prepara- 
tion devices. Thus, the remarkable POVM elements of 
the PNRHD can be used not only to optimally detect 
the appropriate states [7] but also to prepare them from 
quadrature-entangled beams [8j. 
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As an experimental demonstration of the PNRHD 
capabilities, we have tomographically reconstructed a 
group of classical benchmark states. To our knowledge, 
this constitutes the first full tomographic reconstruc- 
tions involving PNR detectors. Indeed, state tomogra- 
phy based on photon-counting detection has been the 
subject of much theoretical work [9 and the few exper- 
imental implementations reported to date involve only 
binary (on/off) detectors [10 . This is due, in part, to the 
relative infancy of PNR detector technology, which is an 
active area of research with several different approaches 
to photon counting [11 . In the experiments presented 
here time- multiplexed PNR detectors are used, but the 
techniques can be readily extended to other PNR detec- 
tors. The ability to change the measurement basis of the 
detector enables applications other than state reconstruc- 
tion for this detection scheme, such as in non-local state 
preparation, precision quantum metrology [12 , or the im- 
plementation of a continuous- variable entanglement wit- 
ness [13 . 

The POVM elements are derived from an analytical 
model by first considering ideal PNR detectors, able to 
resolve n photons (see Fig. [l](a)). In this case the prob- 
ability of obtaining measurement outcome f3 = {ua^nb) 
for LO setting 7 is related to p by jl4j 

Pf3^ = TTab[U^abU^na){na\a \nb){nb\b], (1) 

with U = e*^(^^a+a^S) unitary operator representing 
the BS, R = cos^(^) the LO coupling, aab = \(^){a\a^ Pb 
the two-mode input state and \na(b)) the photon number 
states to be detected at l^a(6)- Using cyclic properties of 
the trace, Eq. ^ can be written asp/3^ = TT^lpbllpj]. For 
the case of ideal detectors the POVM element IT^^^^ is a 

projector flf^''^ = \x){x\b^ where \x)b = {(^a\U^na)\nb)b. 
For R = 1/2 this can be expressed as 



\x)b = e - ^ / / ^ \ . . — |Q)^- 



(2) 



In any realistic scheme one has to account for unavoid- 
able imperfections in the PNR detectors such as loss, 
non-unit efficiency, and overlap in detector responses. For 
a given PNR detector, this can be done by considering 
the detector design [4], or by experimental detector to- 
mography [15]. The time multiplexed detectors (TMDs) 
used in our experiments accept states of light contained 
in pulsed wave-packet modes. Each incoming pulse is 
split into several spatial and temporal modes by a fiber 
beam splitter network that are subsequently registered by 
avalanche photo diodes (APDs). Because the TMD alone 
is not phase sensitive, its operation can be described as 
a map from the incoming photon-number distribution p 
(the diagonal components of the density matrix) to the 
measured click statistics k hy k = CLp. Here L and C 
are matrices accounting for loss and the intrinsic detec- 
tor structure [4 , respectively. To calculate the POVM 




FIG. 2: Experimental set up for the photon number resolving 
homodyne detector (PNRHD). The output of a Ti:Sapph laser 
is split into two arms (by a PBS), which eventually correspond 
to signal p and LO beams. A piezo moves a mirror to set the 
phase of the LO. Both arms are then recombined and sent 
through the configurable homodyne detector. 



elements implemented by our configurable PNR homo- 
dyne detector, the POVMs for TMD detectors Da and 
Db are determined from their CL matrices (characterized 
by independent methods). The n-th element of the a(b) 
TMD POVM then replaces the projectors \na(b)){i^a{b)\ 
in Eq. ([T]) to obtain the final expression for the POVM 
elements II/37. We note that our TMDs can resolve up 
to 8 photons, setting the number of possible outcomes 
to 81 and truncating the operator Hilbert space to 9 x 9 
matrices. 

By adjusting the local oscillator amplitude and phase, 
the detector POVM elements tlp^ can be tuned to 
project onto different bases. We illustrate this by cal- 
culating 11/3^ for two particular examples, shown in Fig. 
[1] For a sufficiently large LO amplitude and detector effi- 
ciency (7^), the joint click outcomes of the PNR detectors 
can nearly project onto a single-mode Schrodinger-kitten 
state or quadrature-squeezed state, a feature which sug- 
gests our detector can be used both for preparation and 
direct detection of such quantum states. Figure [l] depicts 
the Wigner representation of two normalized POVM ele- 
ments H^™ for a = 1.5, r]a,b = 90% [16, and 50:50 beam 
splitter ratio, where (x,p) label the phase-space conju- 
gate variables. Note that a normalized POVM element 
is a positive definite operator with unit trace, allowing 
it to have a phase-space representation similar to that of 
a density matrix. To quantify the state preparation (de- 
tection) efficiency for these detector configurations, we 
calculate the overlap between the target state p^^^ (i.e., a 
Schrodinger-cat-like state or squeezed state) and the nor- 
malized POVM projector H^™, hy p = Tr[n^™p*^^]. In 
our numerical simulations, p ^ 80% can increase to val- 
ues of up to 90% for unit efficiency PNR detectors. 

To demonstrate the ability of a PNRHD we experimen- 
tally reconstruct various states derived from a coherent- 
state laser pulse using TMDs for the PNR detectors. The 
experimental set up, shown in Fig. [2j consists of three 
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main components: the state and local oscillator prepara- 
tion, and the configurable PNRHD. The signal and local 
oscillator are formed by splitting 90 fs TiiSapph laser 
pulses, centered at 784 nm and cavity dumped to a repe- 
tition rate of 250 kHz. In the time- multiplexed detection 
scheme the APD dead-time (t^ ^ 50 ns) in conjunction 
with the total number of temporal modes infiuences the 
maximal detection rate. We chose a time delay between 
temporal modes of 100 ns after which the after-pulsing 
probability of the APDs drops below 0.1%. This, in com- 
bination with electronic time gating, makes after-pulsing 
effects negligible and reduces dark counts to less than 
5 cnt/sec. The input laser mode is split equally into two 
optical paths by a half-wave plate (HWP) and polarizing 
beam splitter (PBS), corresponding to V and H polar- 
izations, respectively. The LO path has a HWP and PBS 
to independently control its amplitude, which is typically 
set to a small percentage of the total signal. A pair of mir- 
rors in the LO path, placed on a piezo-electric-controlled 
translation stage, set the relative phase {0) between the 
LO and signal. The two beams are recombined into a sin- 
gle path at a second PBS. A beam sampler (BS) with low 
refiectivity followed by a set of calibrated neutral density 
(ND) filters completes the state preparation and directs 
the signal and LO to a HWP-PBS combination that acts 
as the homodyne BS in Fig. [l] (a). The outputs of the 
final PBS are sent to TMDs (Da^b)- The joint detection 
events are collected using a field-programmable gate ar- 
ray (FPGA) interfaced with a computer. The joint statis- 
tics are monitored for 100 different phase settings with a 
total measurement time of 15 minutes. The relative phase 
can be monitored using the light transmitted through 
the BS, which is then interfered using a HWP and PBS 
and detected at a photo diode (PD). The typical fringe 
visibility measured at this phase- monitoring step, which 
indicates the mode match between the signal and LO, 
was approximately 70%. This relatively low contrast is 
due to chromatic dispersion in the beam paths and can 
be improved by placing narrow-band interference filters 
in the initial laser beam path. Material imperfections 
and inaccurate calibration in the homodyne HWP-PBS 
system could in principle be accounted for by the model 
parameter i?, however such deviations were found to be 
below 1% in our optical elements and are not considered 
a significant source of experimental uncertainty in the 
error analysis. 

The local oscillator amplitude |a| = (nLo)^^^ is ob- 
tained by measuring the average photon number in the 
LO beam (^n^^^^) and then multiplying by the transmis- 
sion (T) of a set of calibrated neutral density (ND) filters, 
so that (T(n"^^^"))i/2 ^ (nLo)^/^ with relative error of 
approximately 5%, due to the inaccuracy in T. The detec- 
tor efficiencies {r]a{b)) are obtained by fitting the param- 
eters 7^a(5) in the loss matrix L which retrieve an average 
photon number (^^(6)) ^ (nLo)/2, measured at detector 
I)ci(6) when the signal is blocked. For our experiments 
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FIG. 3: Wigner function and corresponding contour plot (in- 
set) of the experimentally reconstructed phase-averaged weak 
coherent states, with average photon numbers (n^^*) equal 
to (a) 0.29 zb 0.03, (b) 0.61 ± 0.06, and (c) 0.71 ± 0.07; (d) 
one extreme POVM element used in the error estimation 
n™^^, (e) difference between two extreme POVM elements 

ndifT fVmax fYmin 

the fitted detector efficiencies were rja = 0.10 ± 0.01 and 
r]b = 0.15 ± 0.02, where the errors {sr^ ^ 10%) are ob- 
tained by propagating the uncertainty in (nLo)- The low 
efficiency, as compared to the APD quantum efficiency 
(^ 60%), is due to the single- mode fiber network used 
to implement the PNR detection. To estimate the effect 
of these errors in the tomographic reconstructions, we 
built two POVMs {n™"'^^™'")} using the maximum (min- 
imum) possible values around the mean for the set of 

, // max(min)\ max(min) max(min)\ -ttj n i 

parameters ((^lq /^Va \% )• We round 

that, while a change of up to 10% in (riLo) does not 
affect the reconstruction (due to the relatively small am- 
plitude of the LO, as compared to the signal), a change 
of 10% in the efficiencies propagates into a relative error 
in the final estimated average photon number (n®^*) of 
~ 10%. In order to estimate the most-likely state (p®®*) 
that is compatible with the empirical photo-count statis- 
tics p^^^, we use a recursive least-squares algorithm to 

minimize (over p^«*) Y.^^iv^'' - Tr[n;37P'''])', subject 
to the constraints p^"* > and Tr[p^"*] = 1 [7 . 

The first class of states to be examined were the phase- 
averaged coherent states. Here 20 evenly-distributed LO 
phases 9 are chosen between and 27r for the data ac- 
quisition. Figure [3] (a)-(c) show the experimentally re- 
constructed Wigner functions and corresponding con- 
tour plots for three different phase- averaged coherent 
states. The Wigner functions are rotationaly symmet- 
ric and centered about the origin, as expected for such 
states. To quantify the error in the estimated state we 
calculate the variance A between the two extreme es- 
timated states pmax(min)^ obtained using the two ex- 
treme POVMs {n^^^^^™^}, by A = 1 - F, with F = 

|Tr(^(V^(^^^p"^^''v^p™)))P the fidelity between the 
two extreme estimated density matrices. Figure [3] shows 
(d) an extreme POVM element for outcome /3 = (1,0) 
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FIG. 4: Wigner function (left) and amplitude of 
the density matrix in the photon- number basis \pnm\ 
(right) for experimentally-reconstructed weak coherent 
states. Here (x,p) labels the quadratures and {n,m) labels 
the photon- numbers. Insets show corresponding contour plot 
(left) and diagonal matrix elements (right). Rows correspond 
to an average photon number (n^^*) of (a) 0.29 ± 0.03, (b) 
0.59 ±0.06, and (c) 1.53 ±0.1. 

and (e) difference between ft^^^^ - A^^". Such POVM 
elements do not look exactly like a projector onto a 
single photon state. This due to the low detector ef- 
ficiency which mixes POVM elements corresponding to 
higher photon numbers. Additionally, as a second fig- 
ure of merit, we quantify the proximity of the recon- 
structed state p®®* with an ideal phase- averaged coherent 
state p^"^ of average photon number (n®®*), using the fi- 
delity F, by replacing pmax(min) ^PA(est) r^^^ variances 
and fidelities (A,F) resulted in (a) (0.004,0.985), (b) 
(0.010,0.973) and (c) (0.027,0.954), respectively. Note 
that the error increases with the average photon num- 
ber, an expected effect which is due to the truncation of 
the Fock state space present in the modeled POVMs. 

Next, in order to test the phase sensitivity of this 
detector, several coherent states were tomographically 
reconstructed. Figure |4] shows the experimentally re- 
constructed Wigner function W{x^p) and density-matrix 
amplitude for coherent states of different aver- 
age photon numbers. The variance and fidelity are 
used again as figures of merit of the state reconstruc- 
tion, with p^^ the density matrix for an ideal coher- 
ent state with average photon number (n®®*) and phase 
0^^^ taking the place of p^^, (n®^*) was obtained from 
the diagonal elements of p®®* and 0^^^ was obtained as 
an average over the off-diagonal matrix elements using 
Pn,m = (n^«t)(^+^)/2e-(<^"'))e-^(^-^)^"VvnJ^. We 
found (A,F) equal to (a) (0.002,0.999), (b) (0.010,0.998) 
and (c) (0.010,0.997), respectively. 

In conclusion, we have introduced a highly adaptable 
homodyne detection scheme with PNR detectors. This 
approach bridges the gap between phase-sensitive weak- 



field homodyne techniques and PNR detection, opening 
a new avenue of research and realm of applicability for 
PNR technology. The various states the detector can di- 
rectly project onto (squeezed states, Schrodinger cats, 
and displaced Fock states) make it useful not only for 
state detection, but also for state preparation, a direc- 
tion which is currently being explored at our group. The 
detector has proven to be effective for state tomography 
as demonstrated by the experimental reconstruction of 
weak coherent and phase- averaged weak coherent states. 
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